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OUTFLOW, INFALL AND PROTOSTARS IN THE STAR FORMING CORE W3-SE 
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ABSTRACT 

We report new results on outflow and infall in the star forming cores W3-SE SMA-1 and SMA- 
2 based on analysis of ~2.5" resolution observations of the molecular lines HCN(3-2), IICO+(3-2), 
N2H+(3-2) and CH30H(52, 3-41.3) with the Submillimeter Array. A high- velocity bipolar outflow 
originating from the proto-stellar core SMA-1 was observed in the HCN(3-2) line, with a projected 
outflow axis in position angle 48°. The detection of the outflow is confirmed from other molecular 
lines. An inverse P-Cygni profile in the HCN(3-2) line toward SMA-1 suggests that at least one of the 
double cores accretes matters from the molecular core. A filamentary structure in the molecular gas 
surrounds SMA-1 and SMA-2. Based on the SMA observations, our analysis suggests that the double 
pre-stellar cores SMA-1 and SMA-2 result from fragmentation in the collapsing massive molecular 
core W3-SE, and it is likely that they are forming intermediate to high-mass stars which will be new 
members of a star cluster in the W3-SE region. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics - ISM: molecules - ISM: jets and 
outflows - radio lines: ISM - stars: formation 



1. INTRODUCTION 

One of the key questions in the study of star forma- 
tion is how protostars accrete material from their parent 
molecular clouds. Accretion around protostars with in- 
termediate and high mass is not well understood due 
to lack of high-resolution observations. Observation of 
infall is needed to provide direct evidence for accre- 
tion. Low resolution observations of a blue-dominated 
double-peaked spectral profile as expected from a model 
of self-absorption in infalling g a,s clumps (NGC 1333- 
IRAS 2 as an example; see iWard-T hompso n et al 



1996; Ward-T ho^pson fc Bucklevll200lH Jorgense n et al 



2004at IPeretto et al.ll200i r do not unambiguously iden 



tify infall because massive cores usually contain bipolar 
outflows, disk rotation and multiple sub-cores in addition 
to gas infall. Complex dynamical processes can produce 
characteristics similar to those expected from the self- 
absorption model of a molecular core with infall. It is 
crucial to carry out high-angular and spectral resolution 
observations along with comprehensive modeling of the 
dynamical processes involved in a massive star forming 
region. 

There is growing evidence that massive stars form 
in groups, but accretion modeling has been for mate- 
rial c oncentrated in monolithic cores (M cKee fc Ostrikeil 
|2007( ) . The formation of binary and multiple stellar sys- 
tems, especially those with two comparable-mass pro- 
tostars and a large sep aration, is poorly understood 
(|McKee fc Ostrikeil [20071 ). Detailed high-resolution ob- 
servations of a binary/multiple proto-stellar system will 
provide valuable data for understanding accretion and 
outflow from each of the sub-cores in the system. 

W3-SE is a dense molecular core located SE of the 
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well-known massive star formation region W 3-Main. We 
adop t ed a distance o f D = 2 kpc for W3-SE (jBlitz et al.l 
119821; iXu efall 120061 ). The spectral energy distribution 
(SED), along with a. double co ntinuum core revealed in 
our previous paper ()Zhu et al.ll20r0( ) suggests that a to- 
tal mass 65 Mq for the W3-SE core over a scale of 10" 
is mainly composed of two proto-stellar cores of compa- 
rable masses, which are separated by 3.5" or 7000 AU as 
revealed in the continuum observations with the Submil- 
limeter Array (SMA) Observations of HCO+(1-0) using 
the CARMA telescope at 6" resolution show that out- 
flows appear to dominate the observed spectra although 
infall might also contribute to the spectrum toward the 
main molecular core which coincides with the double 
cores. Spectral-line observations with higher angular res- 
olution are necessary to separate the contributions of in- 
fall and outflow from the double cores in W3-SE and to 
investigate the possibility of a protobinary system. 

In this paper, we present results of analysis of spectral- 
line data observed with the SMA at ~2.5" resolution. 
Section 2 describes the observations and data reductions. 
Section 3 shows the results on the basis of our analysis 
of the molecular line data including HCN(3-2), HC0+(3- 
2), N2H+(3-2) and CH30H(52,3-4i,3). In Section 4, we 
present an interpretation of the observations in the con- 
text of outflow and infall. Section 5 gives a summary and 
conclusions. 

2. OBSERVATIONS 

The observations toward W3-SE consist of two tracks 
obtained at A 1.1 mm using the SMA in the compact 
configuration with an angular resolution of ^2.5" in 2008 
October 27th and December 20th. The observations were 
centered at R.A.(J2000) 02'^25'"54^50, DEC(J2000) 
= 62°04'11".0, which is used as the reference position 

* The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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in the following analysis. The LO frequency (t'Lo) for 
the 2008 October track was tuned to 271.7 GHz and i^lo 
= 285.3 GHz for the other. The observations consist of 
a lower side band (LSB) and an upper side band (USB) 
separated by 10 GHz, with a bandwidth of 2 GHz for each 
sideband. The typical system temperatures (Tgys) are 
180-200 K and 200-250 K for the observations in October 
and December, 2008, respectively. 

Data re duction and ima ging; were carried out in 
MIRIAD (|Sault et all I1995D with specific implementa- 
tions for SMA data reductior0. The bandpass was cali- 
brated using QSOs 3C454.3 for the lower- frequency track 
and 3C273 for the higher-frequency track. For both 
tracks the complex gains were calibrated using QSO 
0244-f 624. The flux density scale was bootstrapped from 
Uranus using the SMA planet model. By combining the 
visibility data from both LSB and USB, the continuum 
images achieve rms noises of 1.5 mJy and 3.0 mJy for the 
data observed in 2008 October and December, respec- 
tively. Spectral line images were made by subtracting 
the continuum emission from each visibility data using 
the MIRIAD task UVLIN. 

3. RESULTS 
3.1. HCN(3-2) 

A transition of hydrogen cyanide, HCN(3-2) at fo 
= 265.886 GHz, was observed in the LSB of the 2008 
October observations. Figure la presents the HCN(3- 
2) channel maps, showing multiple velocity components 
distributed in the W3-SE core. In the velocity range 
between -39.3 and -38.0 km s^^ a gap with no sig- 
nificant line emission signals is observed. The figure 
also shows compact HCN(3-2) components in the high- 
velocity channels (< -46.7 km s~^ and > -30.5 km s~^). 
Figure 2a shows the integrated line intensity map of 
HCN(3-2), suggesting the presence of three emission 
peaks from the HCN(3-2) line within the main molec- 
ular emission component A, a.s observed by C ARMA at 
an ang ular resolution of 6" (|Zhu et al.l [20100 . We de- 
note the three sub-components as A-NE, A-SW and A- 
SE, respectively. The HCN(3-2) line components A-NE 
and A-SW are located 2" northeast and 4" southwest of 
the proto-stellar core SMA-1 while the component A-SE 
nearly coincides with SMA-2. The integrated intensity of 
HCN(3-2) line emission in A-NE and A-SW is stronger 
than that from A-SE. 

Figure 3 shows the spectra of the HCN(3-2) line in the 
core region of W3-SE with a boundary outlined by the 
square in Figure 2a. A significant absorption feature is 
shown in the core region at velocities in the range be- 
tween -39.3 and -38.0 km s~^. The component A-NE 
appears to be associated with a broad blue-shifted line 
wing, while A-SW shows a broad red-shifted line wing. 
The highly blue-shifted line wing associated with A-NE 
shows a maximum absolute radial velocity of 35 km s~^ 
(above the 5 a cutoff), with respect to SMA-1. Here an 
LSR systemic velocity V^ys = -39.1 km s~^ is used, which 
is determined from the peak velocity of the optically thin 
fine CH30H(52,3-4i,3)- In the region of A-SW, the red- 
shifted line wing shows a maximum absolute radial ve- 
locity of 33 km s^^ (above 5 a) with respect to SMA-1. 

^ http:/ /www. cfa.harvard.edu/sma/miriad 



The velocity profile for both blue- and red-shifted broad 
line wings can be fitted to a power law (~|F — T^ysl"''') 
with an index of 7 « 2.0. The power-law distribution 
of molecular content as function of velocity observed is 
consi stent with num erical models of bow-shock outfiow 
fe.g.. lDownes fc Rav 1999: L_ee_et_al. 2001). 

In Figure 4a, we show the high- velocity gas in both the 
blue-shifted (-74.1 to -45.2 km s~^) and red-shifted (- 
31.9 to -5.8 km 3"^) outflow lobes, located NE and SW 
of the proto-stellar core SMA-1 with separations of 2.2" 
and 5.0", respectively. The configuration of the emission 
from the high-velocity wings, along with the location of 
SMA-1 suggests a bipolar outfiow originating from SMA- 
1, a possible energy source in the region. The major axis 
of the outflow is shown as a straight line connecting the 
two peaks of the blue- and red-shifted knots in Figure 4a, 
at P. A. = 48°. Although the location of SMA-1 appears 
to be offset from the outfiow axis by ~0.4", about 5% of 
the distance between the two high-velocity knots, such an 
apparent offset could be due to the uncertainty (~0.2") 
in the peak position of SMA-1 given by Gaussian fitting. 

In summary, the high- velocity HCN(3-2) line emission 
from the two compact molecular components A-NE and 
A-SW appears to be excited by a high-velocity bipolar 
outflow originating from the proto-stellar core SMA-1. 

3.2. HCO+(3-2) 

The transition of HCO+(3-2) at 267.558 GHz was also 
included in the LSB of the observations at i^LO = 271.7 
GHz. Figure lb shows the channel maps of the HC0+(3- 
2) line. The overall distribution of HCO+ gas in veloc- 
ity is con sistent with that of HCO+(1-0) observed with 
CARMA (iZhu et al.ll20Tol) and simflar to that of HCN(3- 
2) in the SMA observations, showing multiple emission 
peaks at various velocities and a gap in the velocity range 
between -38.9 and -37.9 km s~^. 

However, not all of the emission peaks observed in 
HCO+(3-2) coincide with those in HCN(3-2). In the in- 
tegrated intensity map of HCO+(3-2) in Figure 2b, the 
emission peak A-NE observed in HCN(3-2) appears to 
be embedded in the surrounding emission, and A-SW 
in HCO+(3-2) is not as prominent as it is in HCN(3- 
2). Instead, strong emission in HCO"'"(3-2) appears to 
be concentrated in a bright ridge 0.5" northeast of the 
line connecting the double proto-stellar cores SMA-1 and 
SMA-2. Outside the main core A, the molecular sub- 
core B and the molecular tail C are also present, in good 
agreement with the image of HCO"''(1-0) observed with 
CARMA. In comparison with the HCN(3-2) image, the 
HCO+(3-2) line appears to trace a more extended enve- 
lope around the double proto-stellar cores. 

As Figure 5 shows, the observed HCO+(3-2) spectra 
are also found with broad line wings at the position of 
SMA-1 and its adjacent region, although the line wings 
in HCO"'"(3-2) are not as significant and broad as that 
in HCN(3-2) spectra. The integrated intensity map for 
the blue- and red-shifted line wings is shown in Figure 
6a. The velocity ranges for the integrations are (-57.9, 
-45.3) km s~^ and (-31.8, -18.8) km s~^, with the same 
lower velocity limits as that for HCN(3-2). In Figure 6a, 
the compact blue- and red-shifted lobes are offset from 
SMA-1 simflar to the case of HCN(3-2), aUhough they 
are closer to SMA-1 (1.3" for the blue lobe, and 3.4" 
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for red lobe) than the lobes in HCN(3-2). This could 
be explained if HCN(3-2) emission is excited close to 
bow shocks, while HCO"'"(3-2) traces colder and more 
extended gas which is farther from bow shocks. An- 
other possible reason is that the emission of the high- 
velocity outflow lobes might be blended by the HC0"*'(3- 
2) emission from a large-scale outflow ()Zhu et al.ll201ClD 
and other dynamical components associated with SMA- 
2. Indeed, in Figure 6a a red-shifted HCO+(3-2) com- 
ponent is found close to SMA-2 at high velocity. The 
same component is also marginally detected in HCN(3- 
2) (see Figure 4a), and it might be an outflow component 
associated with SMA-2. 

The low- velocity line wings of the HCO+(3-2) spec- 
tra were also integrated with the velocity ranges (-44.9, 
-41.7) km s^i and (-36.3, -32.3) km s'^ for blue and 
red lobes, respectively. The velocity ranges have sim- 
ilar lower velocity limits as the o nes used in previous 
HCO+(1-0) study (IZhu et al.ll201(l . As Figure 6b shows, 
the distribution of blue- and red-shifted low-velocity 
HCO'^(3-2) gas is consistent with the results of HC0+(1- 
0), implying a possible large-scale outflow in a NW-SE 
direction. 

3.3. N2H+(3-2) 

A transition of protonated nitrogen, N2H+(3-2) at i/q — 
279.512 GHz was included in the LSB of the 2008 Decem- 
ber observations. Figure 2c shows the integrated inten- 
sity map of N2H+(3-2), indicating muhiple N2H+(3-2) 
emission clumps which form a curved filamentary emis- 
sion region. Among the numerous N2H+(3-2) emission 
peaks, the two strongest emission peaks appear to be lo- 
cated near SMA-1 and SMA-2 with a displacement of 
'^O.S" from the continuum peaks, suggesting that a close 
relationship between the dense molecular gas and the 
proto-stellar cores. 

The bright N2H+(3-2) emission region is located SW 
of the double continuum cores, opposite to the locations 
of the bright emission ridge of HCO+(3-2). Emission 
at the locations of SMA-1 and SMA-2 might be influ- 
enced by excitation or chemical conditions. The primary 
formation routes for both HCO+ and N2H+ involve re- 
actions with H;^. For a higher CO abundance, the main 
formation route is the reaction between CO and in 
forming HCO^. While when the CO abundance gets 
low enough, the reaction between H^j' and N2 in form- 
ing N2 H^ becomes the main rn echanism for the removal 
of (| Jorgensen et al.l r2004bl ) . In addition, as CO re- 
turns to the gas phase from a freezing-out phase, destruc- 
tion of N2H"'" through reactions with CO becomes the 
dominant removal mechanism for N2H+. The possible 
anti-correlation of the bright emission regions between 
N2H+(3-2) and HCO+(3-2) near the proto-stellar cores 
might reflect a gradient of the CO abundance across the 
W3-SE core. 

The line profiles of N2lI+(3-2) spectra towards SMA-1 
and SMA-2 are shown in Figure 7. We noticed the skew- 
ness of the line profiles, i.e., at the position of SMA-1, 
the N2H+(3-2) spectrum mimics a "blue-profile", while 
toward SMA-2 N2H+(3-2) is like a "red-profile". The 
main hyperfine structures of the transition J=3-2 for 
N2H'*" spread over a frequency range of only a few 100 
kHz (100 kHz - 0.1 km s'^), while the velocity differ- 



ence between the two emission peaks of the N2H+(3-2) 
spectrum toward SMA-1 is ^5.0 km s^^. Therefore, it 
is unlikely that the observed skewness in the N2H'*'(3-2) 
line profile is due to hyperfine structure. 

In order to model the overall kinematics of the molec- 
ular filament in the region, a centroid velocity map 
of N2H'''(3-2) was made (see Figure 8a) by calculating 
Ween = Jvl-udv/ J lydv at each position. Here Vccn is 
the centroid velocity (intensity- weighted mean velocity), 
and ly is the intensity of the spectrum at a given veloc- 
ity V. In general, the distribution of the centroid velocity 
shows a gradient of 30±10 km s~^ pc~^ over a spatial 
scale of ~20", i.e., the gas in NW is red-shifted while 
that in SE is blue-shifted. This velocity structure ap- 
pears to be similar to the large-scale velocity structure 
traced b y HCO+(1-0) th at is believed to be a result of 
outflow (|Zhu et al.ll2010l ). 

3.4. CH30H(52,3-4i,3) 
The methanol transition CH30H(52, 3-41,3) at = 
266.838 GHz was observed in the LSB in the October 
2008 observations. Figure 2d shows the integrated inten- 
sity map of this transition. The two main CH30H(52,3- 
4i^3) emission peaks coincide with the two proto-stellar 
cores SMA-1 and SMA-2 within the uncertainty in po- 
sition (~0.2"), indicating that this line traces the dens- 
est part of the cores and could be used to determine 
the systemic velocity. The methanol emission appears to 
mainly arise from the two proto-stellar cores as well as 
the high- velocity outflow from SMA-1. In addition, in 
region B, a weak signal of CH30H(52,3-4i_3) is detected 
at a 3cr level. No significant emission from CH30H(52,3- 
4i,3) was detected towards region C. The distribution of 
the centroid velocity in CH3OH emission (Figure 8b) has 
no significant velocity gradients, although at the north- 
west boundary and the position of A-SW the centroid 
velocity is slightly red-shifted. 

3.5. Comparison of Molecular Line Spectra towards 
SMA-1 and SMA-2 

Figure 7 shows the spectra from the peak positions of 
SMA-1 and SMA-2 in the lines of HCN(3-2), HC0+(3- 
2), N2H+(3-2) and CH30H(52,3-4i,3), respectively. In 
general, double-peak profiles appear to characterize the 
optically thick lines toward both cores from different 
molecules. For a given molecular line, the double-peaked 
profile of SMA-1 is dominated by the blue-shifted fea- 
ture, whereas in SMA-2, the red-shifted peak is stronger 
than the blue-shifted one. For a given core, the sepa- 
ration between the peaks of the double-peak line profile 
appears to vary. For example, in the case of SMA-1, the 
velocity separations between the two spectral peaks vary 
significantly from 7.6 km 8"^ for the HCN(3-2) line to 4.0 
km s-i for both the HCO+(3-2) and N2H+(3-2) lines, 
suggesting that different molecular lines may trace dif- 
ferent dynamic processes or chemistry in a proto-stellar 
core. For the HCN(3-2) hne from SMA-1, the large ve- 
locity separation of 7.6 km s~^ between the double emis- 
sion peaks is probably due to the high-velocity bipolar 
outflow because the HCN(3-2) line is the best tracer of 
the shocked and compressed medium among the four ob- 
served molecular lines shown in this paper. 

In contrast to the double spectral peaks observed 
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in HCN(3-2) and HCO+(3-2), the hot molecular line 
CH30H(52,3-4i^3) shows a single-peaked profile toward 
both the proto-stellar cores. Gaussian fits to the spec- 
tra show that CH30H(52,3-4i_3) are peaked at -39.1±0.1 
km s"^ for both SMA-1 and SMA-2. We use this as the 
systemic velocity in this paper, i.e. 14ys = -39.1 km s~^. 
The apparent peak velocities of the spectra have a dif- 
ference of ~0.4 km s~^, equivalent to the width of one 
channel. Figure 8b also suggests a difference ^^0.5 km 
between the centroid velocities of CH30H(52,3-4i^3) 
toward the two cores. However, this inconsistency could 
be due to the slightly asymmetric CH30H(52,3-4i_3) line 
profile toward SMA-2 and contamination from the com- 
plex kinematics in W3-SE region (e.g. the outfiow from 
SMA-1). Thus, we conclude that there is no significant 
difference between the radial velocities of CH3 0H(52 3- 
4i,3) from SMA-1 and SMA-2. 

An emission bump at -40.2 km s~^ is present in 
the spectral dip between the double HCN(3-2) emis- 
sion peaks toward SMA-1 (see left-top panel in Figure 
7) . The velocity of the emission bump appears to corre- 
spond to the velocity of the blue-shifted peak in both the 
HCO+(3-2) and N2H+(3-2) hues, which is blue-shifted 
with respect to the systemic velocity. In addition, the 
red-shifted dip between the two main emission peaks ap- 
pears to be an absorption feature and the intensity be- 
comes negative (~ 5 a) at -37.6 km s^^. We will discuss 
these spectral features in Section 4.2.1. 

4. DISCUSSION 
4.1. High-velocity Bow-shock Outflow in HCN(3-2) 

Two possible mechanisms for accelerating molecular 
gas in outfiows from proto-stellar cores are st ellar winds 
and highly coUimatcd jets (e.g. iDownes & Ray 1999). 
The broad velocity distribution with a power law in both 
blue- and red-shifted HCN(3-2) line wings and the mor- 
phology of the high-velocity HCN(3-2) lobes suggests 
that high-velocity HCN(3-2) emission is excited at jet- 
driven bow shocks. An outfiow model with jet-driven 
bow shocks h as been studied both analytically and nu- 
merically (e.g. iDownes fc R,avl[l999l : ILee et al.l ' l2nnnh . In 
the case of W3-SE, a P-V diagram along the major axis 
of the bipolar outfiow observed in the HCN(3-2) line 
(Figure 9) shows similar characteristics to that predicted 
from the bow-shock model. The broad blue- and red- 
shifted velocity gas components show significant offsets 
of 2.2±0.2" and 5.0±0.2" from SMA-1, providing evi- 
dence that the high velocity HCN(3-2) gas is indeed ex- 
cited and accelerated at the bow-shocked sites outside of 
the proto-stellar core SMA-1. 

To characterize the observed high-velocity outfiow, 
we adopted a simp le model of bow-shock outflow 
(jPownes fc Ravl[l999f ). in which the geometry and ve- 
locity of the molecular gas along the bow-shock surface 
can be worked out in analytical forms. Using the model 
illustrated in Figure 10, a cylindrical symmetric outflow 
radius (r) and major axis (z) plane are used following 
IDownes fc Ravi ()1999l ) . We derived equations for the four 
boundaries that outline the observed P-V diagram (see 
Appendix). A power-law (z ~ |r|*) with an index s = 2 
for the locus of the bow-shock surface is assumed in fit- 
ting the model to the observed P-V diagram of the bipo- 
lar outflow from SMA-1 as discussed below. The sys- 



temic velocity, Vsys = -39.1 km s~^, and the projected 
angular offsets of the blue- and red-shifted bow-shock 
tips [ob ~ 2.2" and an sa 5.0") from the outfiow ori- 
gin have been determined directly from the peak veloc- 
ity of the hot molecular line CIl30II(52.3-4i_3) (Figure 
7) and the IICN(3-2) outflow map (Figure 4a), respec- 
tively. The remaining two free parameters a (inclination 
angle) and v (the velocity of bow-shock with respect to 
the origin of the outfiow) were assumed to have equal 
values for the blue- and red-shifted lobes, and they can 
be determined by comparing the observed P-V diagram 
and the curve from the bow-shocked outflow model. Fig- 
ure 9 shows the observed P-V diagram in contours and 
the curves derived from the bow-shock outfiow model in 
thick lines. A shock velocity w = 40 km s""'^ and three 
different inclination angles a = 55°, 70°, and 85° were 
used to generate the curves. It appears that the model 
with V = 40±10 km s'^ and a = 70±5° best fits the 
overall line wings of the P-V diagram. To estimate the 
uncertainty (1 u) in the values of v and a, we added 
small offset to cause the derived model curve to offset 
by one contour interval of the P-V diagram (3 a). The 
inferred shock velocity of the bipolar outflow implies a 
dynamic age of 3x10^ yr for the the proto-stellar core 
SMA-1, which is an order of magnitude younger than 
that inferred for t he possible large -scale molecular out- 
flow in the region (iZhu et al.ll2010[ ). indicating that the 
source(s) embedded in SMA-1 is probably very young. 
On the other hand, SMA-1 is know n to be associated 
with mid-IR source (jZhu et al.ll2010t ) which implies the 
existence of heating source in SMA-1. We offer two possi- 
ble explanations, i.e. the protostar in SMA-1 has a large 
mass and evolves more rapid than classic low-mass stars, 
or there are multiple YSOs embedded in SMA-1 and the 
compact outflow and the mid-IR source are actually as- 
sociated with different YSOs. 

4.2. Infalling Gas toward SMA-1 
4.2.1. Possible Infall Signature 

The blue-dominated double-peaked proflle of the 
IICO+(1-0) spectrum toward core A observed using 
CARMA at 6" resolution was fltted u sing a model wit h 
both outflow and infall components (jZhu et al.l 120101 ) . 
The fit suggests that outflow significantly contributes 
to the spectral proflle, particularly in the high-velocity 
wings. The high-velocity outflow lobes from SMA-1 are 
separated from the proto-stellar core in the SMA IICN(3- 
2) observations at (2.5" x 2.3") resolution. Contamina- 
tion from the high- velocity outflow in the spectra toward 
the cores has been reduced to a level where an inverse 
P-Cygni proflle, a signature of infall, would become no- 
ticeable if it exists. 

The HCN(3-2) spectra (Figure 3), show blue-shifted 
emission at -40.2 km s~^ toward SMA-1 and the adja- 
cent region, with a sharp dip in the emission at -37.6 km 
s~^. The systemic velocity of -39.1 km lies between 
the two spectral features. The complexity of the HCN(3- 
2) line profile toward SMA-1 might be due to a super- 
imposition of an inverse P-Cygni proflle and a broader 
double-peaked proflle from the high- velocity bipolar out- 
flow. Such an explanation is consist ent with the hyp oth- 
esis that gas is infaUing to SMA-1 C Zhu et al.ir2010f) . 

Figure 7 shows that toward SMA-1 the intensity of 
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the HCN(3-2) dip at -37.6 km becomes negative. 
The spectral dip could be partially caused by unsam- 
pled short-spacings. However, in our observations the 
shortest projected baseline is 5.4 kA, corresponding to 
a spatial scale of ^38". Thus for the spectrum toward 
SMA-1 with a beam size (2". 5x2". 3), flux missing should 
not be a severe problem. We can roughly estimate the 
short spacing effect using the previous study of CARMA 
data. For the HCO+(1-0) spectra at the peak position, 
the flux density of the spectral dip was raised to 0.5 Jy 
beam^^ from 0.0 Jy beam"^ after combining with single- 
dish data. Assuming that the missing flux density of 0.5 
Jy beam~^ is evenly distributed in a structure on a scale 
of the CARMA beam (6".3x5".7), and toward SMA-1 
the flux density of HCN(3-2) is the same as HCO+(1-0), 
the missing flux in HCN(3-2) at the velocity of the dip 
is less than 0.1 Jy within the SMA beam. On the other 
hand, the intensity of HCN(3-2) towards SMA-1 is ~-l 
K (~-0.4 Jy beam"^) at the spectral dip. Therefore, in 
the case of W3-SE the missing flux caused by unsam- 
pled short spacings has only a small influence in the core 
region on a scale of the beam size, and the observed neg- 
ative intensity for the HCN(3-2) spectral dip is likely to 
be real. 

If the emission bump is explained as an inverse P-Cygni 
feature, the infalling gas associated with SMA-1 appears 
to be confined to a region ^6" (^12000 AU) with veloc- 
ities in a range of (-40 to -38 km s~^) that is close to 
the systemic velocity of the proto-stellar core. We made 
an integrated intensity map of this emission in a velocity 
range from -40.0 to -39.5 km s~^ (Figure 4b). This im- 
age clearly shows that the HCN(3-2) feature is associated 
with SMA-1, thus SMA-1 is likely to be the gravitational 
center. Weak HCN(3-2) emission to the SE of SMA-1 at 
a similar radial velocity, could blend with the observed 
HCN(3-2) spectra from SMA-1. 

4.2.2. A General Model Including Both Infall and Outflow 

In order to separate the possible inverse P-Cygni (in- 
fall) feature from the outflow, we adopted a character- 
istic model with a continuum core in the middle com- 
prised by two layers of infalling gas and two layers of 
outflowing gas located outside of the infall c ore, which is 
hereafter referred to as a four-layer model ( Mvers et al.l 
119961: iDi Fr ancesco et aLll200lHAttard et alll2009[) . The 
schematic c onfiguration of the m odel is illustrated in the 
Figure 4 of lAttard et"all (|2009D . 

In the following discussion, the subscripts /, r, B, R 
and c are used to represent the front infall layer, rear 
infall layer, blue-shifted outflow layer, red-shifted out- 
flow layer and the continuum core, respectively. On 
the assumption of homogeneous and isothermal layers 
of gas, the line intensity or brightness temperature from 
the "four-layer complex" can be solved from the radia- 
tive transfer equation for given excitation temperatures 
Tf and of the front and rear layers of infall. Tout of 
both the outflow layers, and Tc of the continuum emis- 
sion. The line brightness temperature, after subtracting 
the continuum in which the contribution from the cosmic 
background radiation is negligible, can be expressed as: 



J{T,)f,{l-e-^) 
J{Tf) 



J(Tout) 



(1 - e-^-)e-^^- 

J(T,)/,(1 - e-^) 



J (Tout) 



fl-e- 



(1 



(1) 



1 ) is a radiation temper- 



(1-e- 



where J(T) 

ature of each layer and /c is the beam filling factor for 
the continuum emission. Equation 1 shows that, in the 
four-layer model, the resultant line-profile can be ex- 
pressed as two separated terms. The term in the first 
bracket on the right side of Equation 1 describes the 
line profile of the in-falling gas, corresponding to an in- 
verse P-Cygni profile in the case of J{Tc)fcTc>J{Tf). 
The infall line profile is attenuated by the absorption 
in the blue-shifted outflow (e""^^). The term in the 
second bracket corresponds to the line profile from the 
bipolar outflow showing a classical P-Cygni profile if 
J{Tc)fcTc>J{Tovit) or a red-dominated double-peaked 
profile if J(Tc)/cTc< J(Tout)/out- The emission from the 
red-shifted outflow layer is attenuated by the absorbing 
gas in both the infall layers and the blue-shifted outflow 
as well as the continuum dust [e^^''^^''^^'^^'^"). 

Each of the optical depths in the four layers can be 
modeled using a Gaussian function with a mean velocity 
and a velocity dispersion for infall (V^n, fin) and out- 
flow {Vont, Cout) assuming the systemic velocity of Vsys- 
Therefore, a total of fourteen parameters, namely, the 
five kinematic parameters {Vi-n, Vont, txin, tJout and T4ys) 
along with four excitation temperatures {Tf, T^, Tout and 
Tc) and five optical depths (ry, r^, tb, and fcTc) are 
required to characterize the physical conditions of the 
gas in each layer. 

4.2.3. Fitting the HCN(3-2) Spectra 

We made a tentative fitting to the observed IICN(3-2) 
spectra with the model described by Equation 1. In the 
case of W3-SE, the quantities for the continuum emission 
[Tc = 40 K and Tc = 0.5) can be de termined from the 
SED as discussed bv lZhu et al.l (j2010( ) assuming that the 
ratio of 2:1 in dust content between SMA-1 and SMA-2 
is proportional to the ratio of the total between the two 
sources. Therefore, thirteen free parameters are needed. 

In order to fit the four-layer model to the observed 
IICN(3-2) line profiles, we first assumed the two out- 
flow layers are optically thin at low velocities and fit the 
low-velocity part (from -43.5 to -34.1 km s""'^), or the as- 
sumed infall part, with the least-square approach. The 
high-velocity parts of the blue wings (< -43.5 km s"^) 
and red (> -34.1 km s^^) of HCN(3-2) are related to 
the bow-shock, and could not be fitted with Equation 
1. Therefore, they are fitted separately with a power- 
law profile (oc \V — Vsys\~'')- The maximum intensities 
of the power-law line wings are set to be at -43.5 and 
-34.1 km s~^ for blue and red wings, respectively, and 
the values are obtained from fitting the low-velocity line 
profile with Equation 1. 

The best fitting results are shown in Figure 11a, which 
gives an infall velocity 0.9 km s~^. The derived excitation 
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temperatures are 6~^\ and 11±3 K for the front and rear 
layers, respectively. The best fits for the thirteen free pa- 
rameters in the final stage was adjusted by eye to achieve 
a minimum residual. The errors (Icr) were assessed by 
adding an offset (Scr) to a best-fitted quantity with other 
parameters fixed, to increase the offset of the best-fitted 
curve from the observed profile (or a factor of 3 increase 
in the residual). The higher excitation temperature for 
the rear layer might suggest the presence of a tempera- 
ture gradient in radius (temperature increases as it gets 
close to the protostar). Thus, an observer receives emis- 
sion from colder part in the front and warmer part in rear 
layers, considering the optical-depth effect. On the other 
hand, the excitation temperature of the in-falling gas is 
also constrained by the observed brightness temperature 
in the presence of an inverse P-Cygni profile. As both an- 
gular and spectral resolution improve, observations with 
interferometer arrays at millimeter and sub-millimeter 
wavelengths will place more meaningful constraints on 
the in-falling gas in proto-stellar cores. For the bipolar 
outflow in SMA-1, the gas components in both the blue- 
and red-shifted outflows are optically thin (rout ~0.1) 
with excitation temperatures of ^75 K, which are much 
higher than the beam-averaged peak brightness temper- 
ature (~8 K) of HCN(3-2). As for the high- velocity parts 
of the spectrum, the power-law index 7 for the blue and 
red lobes of the outflow are 2.0 and 2.5, respectively. The 
difference of the two values could be due to a slight differ- 
ence in the directions of the ejected blue- and red-shifted 
high- velocity gas components. 

Fitting the four-layer model to the spectral profile of 
SMA-2 suggests that a bipolar outflow might also be 
present in this proto-stellar core (see Figure lib). In 
comparison, the expected outflow from SMA-2 has a sim- 
ilar optical depth and excitation temperature to those de- 
rived from SMA-1 but it appears at an earlier stage. The 
observed line proflle in the low- velocity range appears to 
be dominated by outflow emission and no obvious infall 
signatures are present toward SMA-2. SMA-2 is different 
from SMA-1 since it has no obviously associated mid-IR 
point sources. The relative extended mid-IR emission 
in the adjacent of SMA-2 might be due to the external 
heating from the formed star in the W3-SE region. 

We compared these fits with HCN(3-2) data with pre- 
vious fits to CARMA HCO+(1-0) data of lower angu- 
lar re solution (^6") using a similar four-layer model 
(jZhu e t al. 2010). All the parameters from the IICN(3- 
2) fits (see the caption of Figure 11) are consistent with 
the previous HCO+(1-0) fits, especially the infall veloc- 
ity (0.9 km s-i). Different from the HCO+(1-0) flts, the 
IICN(3-2) suggests an inverse P-Cygni profile instead of 
a blue profile, which could be due to the larger contribu- 
tion from the high-velocity outflow in HCN(3-2) and the 
signiflcant improvement of angular resolution. 

We found it is implausible to obtain consistent flts 
with Equation 1 for the entire W3-SE region, and the 
fltted values of t/, t^, Tj and Tr change dramatically at 
difference positions. It could partially due to the con- 
taminations from other components, such as potential 
further multiplicity, the possible large-scale outflow, and 
the possible high-velocity outflow associated with SMA- 
2. A more essential cause could be the over-simplifled 
assumptions of the model itself, which assumes complete 



overlaps of the layers along the light of sight. A more 
detailed model should consider the two dimensional ge- 
ometry of the W3-SE region and include the influences 
of the core SMA-2, the additional HCN(3-2) blue-shifted 
component and other possible components. Although 
due to the limit of the model we can only qualitatively 
analyze most of the fltted physical parameters, an ex- 
ception is the infall velocity (V^n)- It can be consistently 
fltted for the spectra from different positions; its value is 
always between 0.8 and 0.9 km s~^ and it is not sensitive 
to other fltted parameters. Actually it is also consistent 
with the velocity difference between the HCN(3-2) emis- 
sion "bump" and the assumed systemic velocity. 

Although fltting the HCN(3-2) spectrum towards 
SMA-1 is a tentative exercise, the fltting demonstrates 
a possible explanation: a complex including both high- 
velocity wings and an inverse P-Cygini profile, for the 
observed HCN(3-2) line profile toward SMA-1 region. A 
multiple-outflow model for the observed profile toward 
SMA-1 can be ruled out since no absorptions against the 
continuum are expected from the outflow gas with an 
excitation temperature of ~ 75 K, which is an order of 
magnitude larger than the brightness temperature of the 
continuum emission (8 K) at the observing wavelength. 
However, an alternative model with multiple outflows to- 
gether with a cold envelope might fit the observed spec- 
tral profile. High-resolution observations of dense molec- 
ular lines at higher frequencies and a more sophisticated 
model are needed to unambiguously determine the na- 
ture of the observed profile towards SMA-1. 

4.3. The Nature of the Double Cores in W3-SE 

The W3-SE dust core is associated with a total mass 
of 65±10 M0 in a 15" region, and is characterized by 
two continuum emission sources SMA-1 and SMA-2 with 
a projected angular separation 3.5" or 7000 AU. SMA- 
1 and SMA-2 contribute 1/2 and 1/4 of the total flux 
density in 1.1-mm continuum, respectively ()Zhu et al.l 
120101 ). Assuming a uniform distribution for both dust 
temperature and dust-to-gas ratio over the entire W3- 
SE region, a mass of 32±5 M0 is inferred for SMA-1 
and 16±2 Mq for SMA-2. In the vicinity of the double 
proto-stellar cores, at least three additional mid-infrared 
sources have been revealed by Spitzer, indicating that the 
double cores located within or near SMA-1 and SMA-2 
might be a young cluster (jZhu et al.ll20Tol ). Both the 
SMA and CARMA observations of the molecular lines 
reveal a common fllamentary envelope of molecular gas 
surrounding the two main proto-stellar cores, suggest- 
ing that SMA-1 and SMA-2 result from fragmentation of 
molecular gas in the same parent cloud. 

The dust temperature of the W3-SE core region (cov- 
ering SMA-1 and SMA-2) --40 K, indicating that the gas 
in the core region is likely to be heated by both internal 
and external energy sources. Indeed, a group of mid-IR 
sources are located east of SMA-1 and SMA-2, and mid - 
IR emission is also detected in SMA-1 ()Zhu et al.ll2010D . 
The high temperature implies a relatively large Jeans 
mass. Assuming a mass of 65 Mq evenly distributed in a 
15" region (^--^0.15 pc), a Jeans mass could be estimated 
by 
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(|Landll998[ ). where H is the mean molecular weight and 
assumed to be 2.35. The inferred Jeans mass is ~7 Mq, 
far less than the derived masses of SMA-1 and SMA-2. 

The observed infall signature might be due to global in- 
fall and/or inflow toward SMA-1. Although the spectra 
of HCN(3-2) and HCO+(3-2) appear to be dominated by 
outflow, they still present significant red-shifted absorp- 
tion (see Figure 7). The HCO"'"(1-0) spectrum toward 
the W3-SE core, which covers bot h SMA-1 and SMA- 
2, also shows the infall signature (jZhu et al.l 12010 ). A 
global infall in W3-SE is likely to occur as the inferred 
Jeans mass is far less than the total mass of the W3-SE 
molecular core. However, both the spectral map (Figure 
3) and the integrated intensity map (Figure 4b) show 
that the possible inverse P-Cygni profile in HCN(3-2) is 
present toward SMA-1, the major one of the double cores 
in W3-SE. 

SMA-1 and SMA-2 may result from fragmentation of 
a common parent molecular cloud. The mechanisms fo r 
binary formation have been discussed by iTohlind ()2002D . 
The location of SMA-1 and SMA-2 within a common fila- 
ment (Figure 2c) suggests that the two cores result from a 
prompt fragmentation mechanism (Tohlinc 2002). Con- 
sidering the projected separation, d = 7000 AU (3.5") for 
SMA-1 and SMA-2, and their total mass, Mtot w 50 Mg, 
they could be gravitationally bound if their actual sepa- 
ration is comparable to projected one. Assuming SMA-1 
and SMA-2 rotate around their mass center with a dis- 
tance d, the orbital period P = (27r<i^)^/^ (GAftot)"^^^ ~ 
3 X 10"* yr and the relative velocity of the two cores ~2.5 
km s~^. Although no significant difference in the ra- 
dial velocities of the double cores is seen in our observa- 
tions, this could be explained if the orbit is nearly face- 
on as implied by the large inclination angle (^70°) of 
the high-velocity outflow in HCN(3-2) originating from 
SMA-1 (Section 4.1). The W3-SE proto-stellar system 
appears to be very young with an age ~ 3x10^ yr in- 
ferred from the kinematic timescale of the outflow. This 
is an order magnitude smaller than the orbital period 
if the molecular core W3-SE is to form a protobinary 
system, indicating it is still an unrelaxed system. How- 
ever, the possibility that SMA-1 and SMA-2 will relax to 
a gravitationally-bound binary system cannot be ruled 
out. Accretion could decrease the separation of the dou- 
ble cores and help the formation of a bin ary as suggested 
in th e theory of dynamical hardening ()Bonnell fc Batd 
MM)- 

W3-SE can be compared with NGC 1333 IRAS 4, 
which was reported as a binary system consisting of two 
proto-stellar cores 4A (9.2 Mq) and 4B (4.0 M©) with 
an angul ar separation ^30" (^9000 AU) at a distance of 
~300 pc (jSandell et al.lll99"lil . In follow-up observations 
with higher- angular resolution, further multiplicities in 
both the sub-cores 4A and 4B have been revealed (e.g. 
iLav et all [T995t iLoonev et aT1 l2000l[). Given the large 
masses for both SMA-1 and SMA-2 in W3-SE and the 
limited angular resolutions of the existing studies, SMA- 
1 and SMA-2 might turn out to be composed of multi- 
ple sub-cores in higher-res olution observat ions . The ob- 
served adjacent IR sources (jZhu et al.ll2010( ) and multiple 
outflows also suggest further multiplicity in SMA-1 and 
SMA-2. 



Following our previous study with continuum and spec- 
tral lines at low angular resolution, we made high- 
resolution observations of multiple molecular lines to- 
ward the double proto-stellar cores SMA-1 and SMA- 
2 in W3-SE. The two proto-stellar cores are located in 
a molecular filamentary structure which shows an ac- 
tive star-forming region in W3-SE. A high-velocity bipo- 
lar outflow originating from SMA-1 is seen in HCN(3- 
2), with confirmation from other molecular lines, which 
shows a star formation activity in W3-SE. The morphol- 
ogy and the velocity profile of the bipolar outflow fit a 
bow-shock outflow model. SMA-1 and SMA-2 appear to 
result from fragmentation of a collapsing massive molec- 
ular core. Global collapse and gas infall are suggested 
by red-shifted self-absorption observed in both HCN(3- 
2) and HCO+(3-2). An inverse P-Cygni feature observed 
in the HCN(3-2) line toward the proto-stellar core SMA- 
1 suggests that SMA-1 is the major accretion source 
in W3-SE. The double proto-stellar cores in W3-SE are 
ideal targets for further observations with higher angular 
and spectral resolution at (sub-) millimeter wavelengths 
in order to understand the formation of clusters with 
intermediate-mass stars. 

We are grateful to the staff of the SMA and CARMA 
telescopes who make these observations possible. We 
also thank the anonymous referee for his detailed com- 
ments and questions which have significantly improved 
the presentation of these results. We thank Wu Y. for 
the suggestions and discussions in forming the idea of 
this paper, selecting the research targets and providing 
the single-dish data of the targets, suggesting the suit- 
able molecular probes, as well as the discussions in the 
analysis of the results. L. Z. is supported by the SAO 
pre-doctoral program during this research, and he is also 
supported by Grant 10873019 at NSFC when staying at 
Peking University. 
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APPENDIX 

BOW-SHOCK MODEL FOR OUTFLOW 

The bow-shock model (jPownes fc Ravl[T999l) is one of the two popular models that provide mechanisms to accelerate 
molecular gas in an outflow. The observed features of the high- velocity bipolar outflow from SMA-1 appears to show 
the characteristics of a bow-shocked outflow as described in both numerical simulations and analytical models. In 
order to fit the P-V diagram from the observed HCN(3-2) outflow with the bow-shock model, we have worked out 
th e formulas to o u tline the boundary of a P-V diagram from a bow-shock surface based on the analytical model 
of IDownes &: Ravi (|1999( ). In this model the bow-shock discontinuity is described as a rigid surface moving in the 
interstellar medium with a value of bulk velocity v. For instance, in a cylindrical coordinate system, the locus of the 
bow shock in a blue-shifted outflow lobe can be described as: 

z = a- \r l", s>2, (Al) 

where z is an axis along the direction of the bulk velocity of the bow-shock surface, r is radius perpendicular to axis 
z, and s is a power-law index to describe the curvature of the bow-shock surface (see Figure 10). In addition, the 
bow-shocked surface is treated as a rigid body. We made the further assumptions used in the following reduction. In 
the bow-shock rest frame, ambient medium moves towards the bow-shock with a value of velocity v. As soon as the 
ambient molecular gas hits the surface of the bow-shock, the motion of the ambient gas alters the direction that is 
tangential to the bow-shock surface with a value of velocity vi : 



where R : 



Vr 



vi= " -.\I— + R\ (A2) 
dz 

— — s(a — z)'-^^^-'/'*, Vz and Vr are the components of the velocity v projected on the z and r 
dr 



directions (in the bow-shock fixed frame), and a is the projected distance between bow-shocked tip and the origin of 
the outfiow. 

Considering the Doppler shift and the geometry described above, a complete surface of a bow-shocked bipolar outflow 
lobes can be divided into four separated parts, namely a combination of front and rear layers for each of the blue- 
and red-shifted lobes. Each of the four boundaries corresponds to a curve in a quadrant of the P-V space to outline a 
P-V diagram. Therefore, four separate functions v'^^ w^j, v'^^ and v'^^ are used to describe the boundaries of bipolar 
outflow in the P-V diagram. Hereafter, the subscripts B or R represent for blue or red velocities and f or r for front 
or rear layers of a given bipolar outflow. 

Giving a front layer of a blue-shifted bow shock, the observed radial velocity v'^j in the LSR frame at a point of the 
bow shock surface is a function the tangential velocity vi. 



Wgf = i=(— cos a + Rb sina) — wsina + V^jys, (A-3) 

V 1 + ^B 



where V^ys is the systemic velocity of the source originating the bow-shock outflow, and a is the inclination angle, the 
angle between the axis of the outflow and the plane of sky. 

Therefore, substituting v\ with Equation A2, the observed radial velocity at a given point (Equation A3) can be 
re- written as: 



'"Bi = 1 , r>2 {- COS a + Rb sin a) \ — + Rg - w sin a -t- V^ys , 



, V 1 

"Rf ^ -, , n2 (~ COS a - Rr sina)\ — + + w sin a -I- Vsys, 
i + rij^ V to 



V 1 

v'bi- = 1 ^ j>2 ('^"^ " + i?BsinQ!)y - "ysina + Vsys, 



■^Rr = -, ,^n2 i'^^^ a - Rr sin a)J^ + R\ + vsma + Vsys, (A4) 
i -|- it^ V io 

where Rb — ssiaB — z)^^""^^^'^^ and Rr = SR{aR + z)^^"^^^/^'^ . The values of the bow-shock velocities for the blue- 
and red-shifted outflow lobes have been assumed to be the same in Equation A4. This assumption is valid since the 
line widths of the blue- and red-shifted wings are very similar as discussed in Section 3.1. We can further simplify 
the model for the velocity proflles by assuming sr = sb = 2. In addition, the projected distances between the bow 
shock tips and the origin of the outflow (a^ and ur) can be directly determined from the outflow image. The systemic 
velocity Vsys can be determined from the peak velocity of hot molecular lines. For SMA-1, Vsys = -39.1 km s~^ is used 
(see Section 3.4). Thus, the remaining free parameters in Equation A4 are the inclination angle a and the bow-shock 
velocity v. 
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Equation A4 describes the observed radial velocities along the boundaries of a bow-shock outflow as a function of z. 
Given a cylindrical coordinate system with an origin that is the same as that of the outflow (as illustrated in Figure 
10), the position on the sky plane (r', z') can be transferred from the coordinates (r, z) on the bow-shock surface for 
a given inclination angle a: 

z' = z cos a + r sin a 
r' = — zsino; + r cosa (A5) 

Using Equations Al and A5, the sky position z' can be related to the outflow coordinate z for each of the four 
boundaries: 

= ^;cosQ! — (as — z)^^^sinQ! 
z'y^^ = z cos a — {an + z) sin a 

z'^j. = z COS a + (ttB — z) sin a 

Zpj. = z cos a + (a/f + z)"'"^^ sin a (A6) 

Thus, the observed radial velocities (Equation A4) for a given bow-shocked bipolar outflow can be implicitly related 
to the coordinate z' on the sky plane with the Equation A6 by eliminating the argument z. With the four implicit 
functions v' = v'{z'), we can outline a P-V diagram along the projected major axis of a bipolar outflow using the 
bow-shocked model. Note that in the above deduction we worked out the geometrical relation without including optical 
depth eflfect. 
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Fig. 1. — The channel maps of HCN(3-2) (a: left panel) and HCO^(3-2) (b: right panel) with a velocity interval of 1.4 km s ^ . Contours 
are (-1, 1, 2, 4, 6,...)x5(t, a = 0.045 Jy beam"! in both maps with an FWHM beam of 2.5"x2.3" (-65°) as marked at right-bottom. The 
number at left-top in each channel image corresponds to the radial velocity. The green crosses mark the positions of the two proto-stellar 
objects, SMA-1 and SMA-2. 
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RA offset (arcsec) RA offset (arcsec) 

Fig. 2.— Integrated intensity maps of (a) HCN(3-2) in the velocity range of (-70.5, -9.8) km s-\ (b) HCO+(3-2) in (-53.0, -21.0) 
km s-i, (c) N2H+(3-2) in (-43.5, -34.5) km s"! and (d) CH30H(52,3-4i,3) in (-44.5, -33.7) km s'^. Contours are (-1, 1,2,3,4.. .)x3(t Jy 
beam~^ km s~^. The rms noises (1 a) are ~0.7, ~0.7, ~0.5 and ~0.2 Jy Beam~^ km s"'^ for the images (a), (b), (c) and (d), respectively. 
The images (a), (b) and (d) were made from the LSB data observed at j^lo = 271.7 GHz with an FWHM beam of 2.5" X 2.3" (-65°). 
The FWHM beam of the image (c) (USB data) is 2.5" X 2.2" (-55°). Crosses denote the positions of the proto-stellar objects SMA-1 and 
SMA-2. Open circles mark the three peak positions of the HCN(3-2) line emission (A-NE, A-SW and A-SE). The positions of the weaker 
emission components B and C observed in HCO+(1-0) with the CARMA are also indicated. 
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Fig. 3. — A grid of spectra for the HCN(3-2) line emission. The covered region is marked as the black rectangle in Figure 2a. The pair 
of numbers in the parenthesis at top-left in each of the spectra gives the offset in units of arcsecond with respect to the reference center. 
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Fig. 4. — (a) The high-velocity outflow map in HCN(3-2), integrated from blue wing in the velocity range (-74.1, -45.2) km s"^ (blue 
contours) and red wing in (-31.9, -5.8) km s"! (red contours). Contours are (-1,1, 2,3,4.. .)x5cr, 1 cr = 0.28 Jy Beam-^ km 8"^ The 
position angle of the outflow axis (the thick straight line) is 48°. (b) The integrated intensity map of the blue-shifted HCN(3-2) emission 
bump in the velocity range of (-39.9, —39.5) km s~^, which could somehow trace the distribution of possible infalling gas. Contours are 
(-1,3,2,3,4.. .)x3(T, 1 cr = 0.1 Jy Beam"! km s'^. For both the figures the FWHM beam is 2.5"x2.3", P.A. = -65°. For both figures the 
background shows the 1.1-mm continuum distribution, and the wedges are in unit of Jy beam~^. 
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Fig. 5. — A grid of spectra for the HCO+{3-2) line emission. The covered region is marked as the black rectangle in Figure 2b. The pair 
of numbers in the parenthesis at top-left in each of the spectra gives the offset in units of arcsecond with respect to the reference center. 
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Fig. 6. — (a) The high-velocity outflow map in HCO"''(3-2), integrated from blue wing in the velocity range (-57.9, -45.3) km s"! (blue 
contours) and red wing in (-31.8, -18.8) km s~^ (red contours). Contours are (-1,1,2,3,4...) x4cr, la = 0.2 Jy Beam~^ km s~^. The thick 
straight line is the same as in Figure 4a. (b) The low-velocity outflow map in HCO+(3-2), integrated from blue wing in (-44.9, —41.7) km 



s (blue contours) and red wing in (—36.3, —32.3) km s 
For both the figures the FWHM beam is 2.5"x2.3", P.A. 
the wedges are in unit of Jy beam~^. 



(red contours). Contours are (-1,1,2,3,4...) X 5(t, Icf = 0.16 Jy Beam ^ km s ^. 
= —65°. For both figures the background is the 1.1-mm continuum image, and 
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Fig. 7.— The spectra of HCN(3-2), HCO+(3-2), N2H+(3-2) and CH30H(52,3-4i,3) (from top to bottom) at SMA-1 (left panels) and 
SMA-2 (right panels). The vertical dash line in the two figures denotes the peak velocities of the CH30H(52,3-4i,3) lines (-39.1 km s~^), 
which is used as the systemic velocity through the paper. 




Fig. 8. — The centroid velocity maps of (a) N2H+(3-2) and (b) CH30H(52,3-4i_3). The contours represent the integrated intensity of the 
line emission with values of (-1,1,2,3,4...) x3cr Jy beam~^ km s"'^ for the given a = 0.5 and 0.2 Jy beam~^ km s"'^ with an FWHM beam 
of 2.5" X 2.3" (-65°) and 2.5" X 2.2" (-55°) in images (a) and (b), respectively. The color wedge scales the radial velocity in the unit of 
km s~^ . 
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Fig. 9. — The position-velocity (P-V) diagram of HCN(3-2) along the outflow axis as indicated with the straight line in Figure 4. The 
origin in the position axis corresponds to the position of SMA-1 projected on the axis. Contours are 0.2, 0.4, 0.6, 0.8, ... Jy beam-l. The 
angular and velocity resolutions are 2.7" and 0.45 km s~^, respectively. The observed P-V diagram is fitted with the bow-shocked outflow 
model as discussed in Appendix assuming the systemic velocity in the LSR frame (Ksys = —39.1 km s~^), the velocity of the bow-shock 
{v = 40 km s~^), and the spectral index of the power law for the bow-shock surface as function of radius (s = 2). The dashed-dot-dashed, 
solid, dashed lines stand for the model curves corresponding to the three different values of inclination angle («=55°, 70° and 85°). The 
curve with o=70° appears to represent the best fitting. Note that bow-shock outflow model shows a good fit to the high-velocity gas. The 
low-velocity gas is probably disturbed by additional processes other than outflow. For example, the emission in the top-left quadrant in 
the velocities of (-40, -45) km s"'^ and the position offsets of (5", 6") are likely from the infalling gas or the gas in the molecular envelope, 
which is not related to the bow-shock outflow. 
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Observer 

Fig. 10. — The sketch illustrates the geometry of the bow-shock outflow model, which is taken from the Figure 6 in lDownes &: Ravi II1999I) . 
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Fig. 11. — (a) The observed HCN(3-2) spectral profile towards SMA-1 (solid black curve) is fitted with the four-layer model (red curve) 
as described by Equation 1, which is discussed in Section 4.2.3 The best fitting gives the optical depths and excitation temperature for the 
infall gas in front and rear layers {t{ = 5.0, = 0.5, J(Tf) = 2 K, J(Tr) = 6 K), and the outflow gas (tb = 0.08, tr = 0.06, J(Tout) = 76 
K), as well as the dust (tc = 0.5, J{Tc)fc = 7 K), respectively. The best fitted velocities and velocity dispersions are (Vin = 0.9 km s~^, 
(Tin = 0.4 km s~^) and (Vout = 3.6 km s^^, (Tout = 16 km s^^) for the infall and the outflow, respectively. The blue dashed and dotted 
lines represent the best fitted spectral components for the infall and the outflow described in Equation 2, respectively. The line wings are 
not fitted with the four-layer model; instead, they are simply fitted using a power-law profile with an index of 2. The black dashed line is 
the residual of subtracting the best fitted model from the data, (b) The observed HCN(3-2) spectral profile towards SMA-2 is fitted to a 
model with outflow alone (Vin = km s~^) as shown in red curve. 



